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ABSTRACT 

The characteristics of propellant mixing near 
the injector have a profound effect on the liquid 
rocket engine performance. However, the flow 
features near the injector of liquid rocket engines 
are extremely complicated, for example 
supercritical-pressure spray, turbulent mixing, and 
chemical reactions are present. Previously, a 
homogeneous spray approach with a real-fluid 
property model was developed to account for the 
compressibility and evaporation effects such that 
thermodynamics properties of a mixture at a wide 
range of pressures and temperatures can be 
properly calculated, including liquid-phase, gas- 
phase, two-phase, and dense fluid regions. The 
developed homogeneous spray model 
demonstrated a good success in simulating uni- 
element shear coaxial injector spray combustion 
flows. However, the real-fluid model suffered a 
computational deficiency when applied to a 
pressure-based computational fluid dynamics 
(CFD) code. The deficiency is caused by the 
pressure and enthalpy being the independent 
variables in the solution procedure of a pressure- 
based code, whereas the real-fluid model utilizes 
density and temperature as independent variables. 
The objective of the present research work is to 
improve the computational efficiency of the real- 
fluid property model in computing thermal 
properties. The proposed approach is called an 
efficient real-fluid model, and the improvement of 
computational efficiency is achieved by using a 
combination of a liquid species and a gaseous 
species to represent a real-fluid species. 
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It is found that the overall computational time 
can be reduced by one order of magnitude by 
avoiding the stiffness of inverting from the 
dependent variables (pressure and enthalpy) to 
independent variables (density and temperature) 
for multiple-species mixtures. The methodology 
of the proposed model and thermodynamics 
properties of some liquid species calculated from 
the model are reported herein. The accuracy and 
computational efficiency of the proposed efficient 
real-fluid model were evaluated by comparing to 
the numerical result of the developed real-fluid 
model tested previously. 


INTRODUCTION 

With the great demands and competitions for 
cheaper and safer launchers access to the space, 
the design of advanced reusable launching 
vehicles (RLV) is required to not only have higher 
payload-to-cost ratio, but also improve the 
reliability and thus extend the life of the reusable 
launching vehicle. The performance and hardware 
integrity of liquid rocket engines is greatly 
influenced by the propellant mixing near the 
injector. Understanding flow characteristics near 
the injector is thus very crucial in the advanced 
liquid rocket engine design. To obtain the 
information of flow phenomena near the injector, 
researchers have to rely on either the experimental 
study or numerical analyses. Conducting hot fire 
test of liquid injectors not only is too expensive 
and costly for all preliminary design concepts, but 
also encounters some technical problems in 
measuring multi-phase flow properties at the 
elevated pressure condition. Hence, most of the 
experimental tests were performed to measure 
chamber pressure and thrust of the engine. For the 
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numerical analysis of the liquid injector, 
traditionally it has been done by using either 1-D 
analytical tools, such as CICM 1 code, or CFD 
codes with ideal gas correlation, such as the TDK 
code 2 . Though the effect of multiple phases was 
accounted for in the 1-D analytical tools, 
developed based on the cold flow data, the effects 
of chemical reactions and spatial non-uniformity 
were not accounted for. Moreover, if the designed 
operating conditions were out of the range of the 
database, then the validity of the empirical 
correction of the analytical code becomes 
questionable. In addition, the environments of the 
injector flow of liquid rocket engines are three 
dimensional because of multiple injectors and the 
vorticities. Traditionally, the operating conditions 
of liquid rocket engine are below the critical 
temperature and well above the critical pressure of 
propellants. Thus, neither the perfect-gas law is 
valid, nor the surface tension is stable enough to 
distinguish the interface between liquid and 
gaseous phases. Though the CFD methodology 
has been mature enough to effectively analyze 
detailed flowfield of combustion, most of CFD 
codes are using the ideal-gas correlation, and thus 
the effect of multiple phases is not accounted for. 
In our previous studies 3 ' 6 , a homogeneous spray 
approach with a real-fluid property model was 
developed to simulate the spray combustion 
problem. The homogeneous spray approach, 
which assumes equilibrium between liquid and gas 
phases (i.e. both liquid and gas are treated as 
continuums, and there are no momentum and heat 
transfer lags), is suitable for near critical and 
super-critical spray and can be computationally 
efficient without tracking liquid droplets because 
the surface tension of liquid phase either does not 
exist or is negligible. The real-fluid property 
model is employed to calculate thermal properties 
with the inclusion of the compressibility effect and 
the heat of evaporation. 

The homogeneous spray approach with the 
real-fluid model was incorporated into a well- 
tested pressure-based CFD code, FDNS 7 ' 9 , as a 
homogeneous spray CFD model to analyze the 
flowfield near the injector as well as in the 
chamber and nozzle of the liquid rocket engines. 
The homogeneous spray CFD model 3 ' 6 
successfully predicted the benchmark test cases of 
the 2 nd International Workshop on Rocket 


Combustion Modeling (IWRCM) 10 ' 12 , and other 
liquid injector configurations 13 ' 14 . However, the 
developed homogeneous spray CFD model was 
handicapped by the computational deficiency, 
especially for the reacting flow. This deficiency is 
primarily in the iterative matrix inversions for 
calculating mixture density and temperature from 
mixture enthalpy and pressure, because the 
pressure and enthalpy are the independent 
variables in the CFD solution algorithm, while for 
the real-fluid mode the density and temperature 
are the independent variables. The computational 
deficiency is worsened by the stiffness of the 
species changes in reacting flows. An efficient 
real-fluid model is thus proposed to relieve this 
constrain by representing a real-fluid species with 
a combination of a gaseous species and a liquid 
species, where the thermal properties of both 
species can be modeled by the perfect-gas type of 
correlations. With this linearization procedure, the 
enthalpy of the liquid species is a function of 
temperature only with a known pressure. The 
matrix inversion process can thus be avoided, and 
the computational efficiency is greatly improved. 
The proposed efficient real-fluid model was 
incorporated into the FDNS code to compute 
thermal properties of local mixtures with the 
homogeneous spray approach. The numerical 
simulation of the RCM test case, previously 
simulated with the real-fluid model, was repeated 
with the proposed model to examine its accuracy 
as well as its computational efficiency. The 
present study is just a proof of concept that the 
proposed model can be applied to simulate multi- 
phase spray combustion flows. Any deficiency 
found from this study will be improved in the 
future. 


NUMERICAL METHODOLOGY 

The characteristics of spray combustion flows 
are multiple phases and chemical reactions. In 
order to properly account for the effect of phase 
change in spray combustion flows, a generalized 
equation of state (EOS) is required in order to 
accurately compute the thermal properties (such as 
density, enthalpy, speed of sound, and etc.) of the 
liquid propellant over a wide range of operating 
conditions. A typical pressure-density correlation 
for a given species can be plotted as shown in 
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Figure 1, which consists of liquid, gas, two-phase, 
and dense fluid regimes. In our previous studies 3 ' 
6 , a real-fluid model was developed and 
successfully predicted the compressibility effect of 
oxygen and nitrogen at an elevated pressure. In 
the developed real-fluid model, thermal and 
caloric equations of state, vapor pressure, heat of 
vaporization, surface tension, and transport 
properties are modeled with the equations of state 
proposed by Hirshfelder, et al 15 ' 16 (we term these 
the HBMS equations of state) and with 
conventional correlations 17 , for the other 
properties. These correlations are explicit in 
density and temperature, and can be expressed as 
follows. 

HBMS thermal equation of state: 
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where p c , T c , and P c are the density, 
temperature, and pressure at the critical condition, 
while p r and T r are the reduced density and 
temperature which are normalized values with 
respect to the value at the critical condition. H and 
H 0 are the real-fluid enthalpy and ideal gas 
enthalpy of a give species. Z c is the 
compressibility for a given species at the critical 
point. By and R are the coefficients of the thermal 
property polynomial and the gas constant for a 
given species, respectively. However, for most of 
the multiple-species flow solvers either pressure 
and enthalpy are the primitive variables (pressure 
based), or the density and enthalpy are the 
primitive variables (density based). The inversion 
process from the dependent variables to the 
independent variables is very computationally 
intensive, especially for the chemically reacting 
flows where a mixture consists of multiple species 


and species concentrations of the mixture can 
change drastically (stiffness problem). The 
efficient real-fluid model is proposed to represent 
a real-fluid species (S) with a mixture of a liquid 
(S,) and a gaseous (S g ) species, and can be 
expressed as 

S = (\-f)S,+fS g 

where / is the fraction of the gaseous species, and 
the split between and Sg is defined by 

P = 0-f)Pi+fP g 

The density (p ) of the real-fluid species is obtained 
from the developed real-fluid model. One of the 
main features of the proposed efficient real-fluid 
model is that the density-temperature-pressure 
correlations of both the liquid and gaseous species 
are similar to the perfect-gas law such that the 
computational efficiency can be improved. Thus, 
the density of liquid (p) and gaseous (p g ) species 
can be explicitly calculated from pressures and 
temperatures. In addition, for the liquid species a 
compressibility factor ( Z c ), which is a function of 
temperature, and a fictitious molecular weight 
( M wl ) are employed to account for the 

compressibility effect. With a known pressure, the 
fictitious molecular weight, the compressibility 
factor, and the liquid density can be calculated as 


M wJ 
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R u , T s , and p sL are the universal gas constant, 

saturation temperature, and the saturation liquid 
density, respectively. The correlation of specific 
enthalpy between real-fluid species ( h ), liquid 
species (A ; ), and gaseous species (h g ) can also be 
expressed as 

h = (\-f)h 1 +fh g 
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The second main feature of the present model is 
that for a given pressure the specific heat (C p l ) 


and the enthalpy ( H ,) of the liquid species are 
polynomial functions of the temperature only 
similar to those of the gaseous species of the CEC 


code 18 , which can be shown as 




n= 1 


K n=l fl 


A„ are polynomial coefficients of the liquid 
species, and can be obtained by curve fitting the 
real-fluid enthalpy for the temperature below 
either the saturation temperature (at subcritical 
pressure) or the critical temperature (at 
supercritical pressure). As a result, the effect of 
the heat of evaporation will be accounted for by 
increasing the fraction of gaseous species as the 
temperature increase. The last main feature of the 
present model is generating a look-up table for the 
speed of sound and the ratio of the specific heat 
based on the developed real-fluid model. The 
need for the speed of sound and the ratio of 
specific heat is to satisfy the correlation between 
the density and the pressure changes, which can be 
shown as 


density and temperature from mixture enthalpy 
and pressure is avoided. 

The proposed efficient real-fluid model was 
incorporated into the homogeneous spray 
methodology in the FDNS 7 ' 9 code to simulate the 
spray combustion problem. The framework of the 
FDNS code is an elliptic, finite difference Navier- 
Stokes flow solver, which employs a predictor 
plus multi-corrector pressure-based solution 
algorithm. Higher order upwind, total variation 
diminishing (TVD), or central difference schemes 
plus adaptive second-order and fourth-order 
dissipation terms are used to approximate the 
convection terms of the transport equations. 
Various matrix solvers, such as vectorized point 
implicit, conjugate gradient, and generalized 
minimal residual 19 (GMRES), are provided in the 
code such that users can select one for a given 
transport equation. Since the FDNS flow solver is 
a structured code, multi-block, multi-zone grid 
topology was employed in the code so that 
problems with complex geometries can be 
analyzed efficiently. Since the flow solver is a 
pressure-based code, a pressure correction (P') 
equation is derived from the continuity equation 
and can be expressed as: 


Ap = - 
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This correlation can be applied to both density- 
based and pressure-based flow solvers, and is valid 
for liquid phase, gas phase, two-phase, and dense 
fluid regions. In order to apply the proposed 
effective real-fluid model to a flow solver, a pre- 
processor was constructed to generate the thermal 
properties of the liquid species for a given 
pressure. In the present model, the temperature is 
an explicit function of the specific enthalpy at a 
given pressure, and the fraction of the gaseous 
species and the density can be directly calculated 
based on the temperature as well. As a result, the 
computational time can be greatly reduced and the 
time marching step size can be increased because 
the matrix inversion for calculating mixture 


pn+i = p n+p' ; p =r / 2 ; r/=-D p VP' 

y / a 

where the superscripts * and n denote the value at 
the intermediate and previous time steps, 
respectively. D p is the inverse of the matrix of the 
coefficients of the convective terms in the finite- 
difference form of the inviscid equations of 
motion. For the homogeneous spray approach, the 
speed of sound ( a ), and the ratio of specific heat 
(y) of a multi-component mixture can be 
calculated from the real-fluid properties of its 
constituents based on the additive volume method. 
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where n , , a, , y , , R, , (C p ), , a, , and M wi are the 

mole fraction, mass fraction, ratio of specific heat, 
gas constant, specific heat, speed of sound and the 
molecular weight of species “i”, respectively. 


NUMERICAL RESULTS 


In order to utilize the proposed effective real- 
fluid model, thermal properties of the liquid 
species as well as the split between gaseous and 
liquid species have to be generated beforehand. A 
pre-processor was constructed to generate those 
data needed by the present model based on the 
real-fluid model. The thermal properties of 
oxygen, hydrogen and RP-1 generated by the pre- 
processor are shown in Figures 2-13, respectively. 
It can be clearly seen that the thermal properties 
generated by the developed pre-processor accounts 
for the compressibility and heat of evaporation 
effects for these three species properly. 

To examine the accuracy of the present model, 
a LOX/GH2 uni-element shear coaxial injector, 
RCM-3 test case (Ref. 1 2), ' was simulated. 
Previously, the homogeneous real-fluid model was 
employed to simulate this test case, and good 
success was reported in Ref. 3 and 5. Thus, the 
accuracy and computational efficiency of the 
present model with the homogeneous spray 
approach can be evaluated comparing to the 
numerical result of the developed real-fluid model. 
As illustrated in Table 1, this is a super-critical 
spray combustion flow test case. Since surface 
tension is either zero or not well defined for super- 
critical conditions, drops either do not exist or are 
very unstable and do not survive very long. 
Therefore, the homogeneous spray model is an 
appropriate and effective approach for simulating 
this test case. The injector configuration and flow 
conditions for the RCM-3 test case are presented 
in Figure 14. A two-zone mesh system (61x39 
and 301x101) was used to model the injector 
section and the combustion chamber. 


As shown in our previous study 345 , the 
numerical results of the equilibrium chemistry are 
almost identical to those of the finite-rate 
chemistiy model because the reaction rate of the 
H 2 /O 2 combustion at the elevated pressure is 
extremely fast. Moreover, the focus of the present 
study is to verify the accuracy of the proposed 
effective real-fluid model in predicting the 
thermodynamics properties over a wide range of 
temperatures. Therefore, it is more 
computationally efficient to utilize the equilibrium 
chemistry to simulate the super-critical O 2 /H 2 
combustion test case in the present study. The 
contours of predicted OH radical concentrations of 
both the present and the real-fluid models are 
shown in Figure 15. It is obvious the numerical 
results of both models are almost identical. In 
Figure 16, the oxygen concentration of the real- 
fluid model compares to the concentration of 
gaseous and liquid oxygen predicted by the 
present model. The sum of liquid and gaseous 
oxygen of the present model agrees very well with 
the oxygen concentration of the real-fluid model. 

The radial temperature profiles predicted at 
several axial stations are shown in Figures 17-19, 
while the axial temperature profiles at several 
radial locations are shown in Figures 20-22. The 
agreement between the present model and the real- 
fluid model is excellent despite some small 
discrepancy at the location very close to the 
injector. This minor discrepancy is caused by the 
approximation that for the entire flow field all the 
property correlations of the present model are 
obtained based on a given pressure (chamber 
pressure in this case), whereas the pressure near 
the injector deviates substantially from the 
chamber pressure. This deficiency can be 
eliminated if the thermal properties are 
interpolated from property correlations of various 
pressures based on the local pressure. The thermal 
property interpolation for the proposed effective 
real-fluid model is being developed as part of this 
on-going research project. Nevertheless, for those 
multi-phase flow problems, where the pressures of 
the entire flow field do not differ too much from 
the value that the thermal properties for the liquid 
species were generated for, the fixed-pressure 
property approach is very efficient and effective. 
In addition to the model accuracy, more 
importantly the overall computational time is 
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reduced by one order of magnitude because the 
numerical simulation of the present model can be 
conducted with larger time marching step size 
(from 0.1 psec to lpsec). This proves that the 
present model is able to relieve the stiffness 
problem by the utilization of two species 
representing a real-fluid species, as well as 
decoupling the process of computing the density 
and temperature. 


CONCLUSIONS 

A property pre-processor based on the present 
model was constructed to generate thermal 
properties for various liquid species. As can be 
seen from the generated thermal properties, the 
effects of compressibility and the heat of 
evaporation are taken into account by the 
developed property pre-processor. The present 
model does meet the goal of improving the 
computational efficiency of the real-fluid model. 
Meanwhile, the numerical result of the present 
model is almost identical to that of the real-fluid 
model. A minor discrepancy between the present 
model and the real-fluid model occurs near the 
injector, which is caused by the pressure deviation 
from the condition, which the thermal properties 
of the liquid species were generated for. This 
deficiency can be overcome by interpolating local 
thermal properties from a pre-generated property 
correlations based on the local pressure. The 
present study is just a proof of concept that the 
proposed model can be applied to simulate multi- 
phase spray combustion flows. The variable- 
pressure version of the efficient real-fluid model is 
being developed as part of this on-going research 
activity. The result of this development will be 
reported in the future. 
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Table 1. RCM-3 Test Case Operating Conditions 


Conditions 

h 2 

02 

Pressure 

6 MPa 

6 MPa 

Mass flow rate 

70 g/s 

100 g/s 

Temperature 

287 K 

85 K 

Density 

5.51 kg/m3 

1177.8 kg/m 3 

Cp 

15100 J/kg/K 

1660.9 J/kg/K 

Velocity 

236 m/s 

4.35 m/s 

Viscosity 

8.67E-6 kg/m/s 

12.34E-4 kg/m/s 
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Typical Thermal Properties Curve 
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Figure 1 . Typical P-p-T correlation for a given 

species Figure 4. Comparison of oxygen speed of sound 

at various temperatures and pressures 



Figure 2. Comparison of oxygen density at Figure 5. Fractions of gaseous oxygen at various 

various temperatures and pressures temperatures and pressures 
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Figure 3. Comparison of oxygen specific enthalpy Figure 6. Comparison of hydrogen density at 

at various temperatures and pressures various temperatures and pressures 
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Figure 7. Comparison of hydrogen specific Figure 10. Comparison of RP-1 density at various 

enthalpy at various temperatures and pressures temperatures and pressures 




Figure 8. Comparison of hydrogen speed of sound 
at various temperatures and pressures 


Figure 11. Comparison of RP-1 specific enthalpy 
at various temperatures and pressures 
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Figure 9. Fraction of gaseous hydrogen at various 
temperatures and pressures 


Figure 12. Comparison of RP-1 speed of sound at 
various temperatures and pressures 
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Figure 13. Fraction of gaseous RP-1 at various 
temperatures and pressures 
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(a) Real-fluid model 



(b) Efficient real-fluid model 

Figure 16. Comparisons of the oxygen 
concentrations near the injector face 
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GH* 


Critical Pressure 5.04 MPa 1.29 MPa 

Critical Temperature 1 54.6 *K 33 *K 
inlet Temperature 85 *K 287 *k 

Mass Flow Rate 0. 1 Hges 0.07 kQti 


Chamber Pressure * 6 MPa 

inlet Turbulent Kinetic Energy - 0.00375 

U^. Propelianr injection speed 


Figure 14. Configuration and flow conditions for 
the RCM-3 test case 


(a) Real-fluid model 




Comparisons of temperature profiles at 
X/D = 10 (RCM-3 test case) 



(b) Efficient real-fluid model 


Figure 15. Comparisons of the OH concentrations 
near the injector face 
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Figure 1 8. Comparisons of temperature profiles at Figure 21 . Comparisons of temperature profiles at 

X/D = 20 (RCM-3 test case) R/D — 2 (RCM-3 test case) 



Radial Distance (mm) Distance from the Injector (mm) 

Figure 19. Comparisons of temperature profiles at Figure 22. Comparisons of temperature profiles 
X/D = 30 (RCM-3 test case) along the chamber wall (RCM-3 test case) 



Distance from the Injector (mm) 

Figure 20. Comparisons of temperature profiles at 
R/D = 0 (RCM-3 test case) 
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